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The  pumpjet  is  considered  to  be  one  of  the  most  promising 
candidate  propulsors  for  high  speed  underwater  vehicles  and, 
as  a  matter  of  fact,  it  has  recently  been  employed  for  MK-48 
torpedoes,  ALWT--Advanced  Light  Weight  Torpedo,  now  called 
MK-50  and  other  underwater  vehicles.  The  pumpjet 
superiority  over  other  propulsion  devices  is  represented  by 
two  major  factors,  i.e.,  high  efficiency  and  quietness. 

The  pumpjet  is  one  of  few  fluid  devices  which  positively 
utilizes  retarded  wake  flow  and  produces  high  propulsive 
efficiency.  This  peculiar  situation  may  be  understood 
readily  by  considering  the  momentum  equatin  applied  to  a 
control  volume  surrounding  an  underwater  vehicle,  fixed  to 
the  inertial  coordinate  system.  In  the  conventional  pro¬ 
peller,  for  example,  the  velocity  of  flow  comng  into  a  pro¬ 
peller  blade  is  approximately  equal  to  the  vehicle  speed 
since  the  propeller  diameter  is  large  enough  to  enjoy  the 
free  stream  flow.  In  order  for  the  propeller  to  generate 
any  effective  thrust,  it  should  accelerate  the  flow,  the 
ejected  flow  speed  being  faster  than  the  incoming  flow.  if 
one  observes  this  situation  from  the  inertial  frame,  the 
ejected  flow  has  a  finite  positive  flow  speed  against  the 
surrounding  environment.  It  means  that  a  certain  amount  of 
the  energy  imparted  on  the  fluid  by  the  thruster  is  dumped 
in  the  surrounding  water.  On  the  other  hand,  the  pumpjet 
receives  the  retarded  flow  velocity,  slower  than  the  free 
stream  velocity.  In  order  to  generate  a  thrust,  again  this 
flow  should  be  accelerated.  However,  if  the  pumpjet  is 
properly  designed,  the  accelerated  flow  velocity  can  nearly 
be  that  of  the  vehicle  speed.  If  one  looks  at  a  similar 
control  volume,  from  the  inertial  f.ame,  the  ejected  flow 
out  of  the  pumpjet  has  almost  no  absolute  velocity  and  thus 
leaves  hardly  any  jet  wake  after  the  vehicle  passed.  There 
exists  much  less  wasted  energy  in  the  flow  field  after  a 
vehicle  with  a  pumpjet  passes.  This  is  the  major  reason  why 
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the  pumpjet  can  produce  such  high  propulsive  efficiency  such 
as  90%  or  higher  if  it  is  properly  designed. 

Quietness  is  a  guaranteed  aspect  with  the  pumpjet,  as  can  be 
seen  from  its  configuration  (see  Figure  1-1);  a  long  shroud 
completely  surrounding  the  rotor  helps  prevent  rotor  noise 
from  emitting  into  the  outside  flow  field.  Furthermore, 
this  "internal"  flow  machine  has  better  resistance  charac¬ 
teristics  against  cavitation,  resulting  in  quieter  shallow 
water  operation  where  propulsors  are  most  susceptible  to 
cavitation. 

However,  in  order  to  achieve  such  a  high  standard  of  perfor¬ 
mance  there  are  many  penalties  to  be  paid  in  reality.  The 
first  such  penalty  naturally  stems  from  the  pumpjet’s  util¬ 
izing  the  velocity-retarded  wake  flow.  A  typical  meridional 
flow  distribution  at  the  inlet  of  pumpjet  rotor  is  shown  in 
Figure  1-2;  the  velocity  at  the  hub  is  only  30%  of  the  free 
stream  velocity  and  rapidly  increases  to  75%  at  the  shroud 
internal  boundary.  This  large  velocity  gradient  in  the 
transverse  direction  is,  of  course,  built  up  by  the  viscous 
boundary  layer  effect  and  is  one  of  the  key  features  causing 
difficulties  in  design,  fabrication  and  eventually  in 
achieving  the  pumpjet  high  performance. 

When  one  designs  an  axial  or  a  near  axial  pump,  it  is  cus¬ 
tomary  to  distribute  the  blade  loading  from  hub  to  tip  in  a 
forced  vortex  or  a  free  vortex  distribution  method,  such  as 
shown  in  Figure  1-3.  Such  distribution  methods  are  impor¬ 
tant  in  obtaining  as  uniform  a  discharge  jet  behind  the 
rotor  as  possible  to  minimize  the  mixing  loss.  However,  a 
serious  problem  arises  in  attempting  to  implement  either 
forced  vortex  or  free  vortex  loading  distribution  against 
the  flow  field  having  a  large  velocity  gradient,  as  shown 
in  Figure  1-2.  Due  to  the  lack  of  enough  meridional  flow 
velocity  near  the  hub,  the  blade  there  should  be  designed 
to  have  extremely  large  incidence  angle  as  well  as  large 
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camber.  It  is  f<.  r  this  reason  that  the  pumpjet  rotor 
designed  to  date  has  a  distorted  profile  shape  from  hub  to 
tip,  see  Figure  1-4.  If  this  were  a  conventional  propeller, 
the  stagger  angle  would  become  smaller  towards  the  hub  and 
the  camber  would  stay  more  or  less  constant.  However,  for 
the  reason  mentioned  above,  the  pumpjet  blade  stagger  angle 
first  becomes  smaller  up  to  the  midspan  area  but  becomes 
larger  toward  the  hub  and  thus  the  camber  is  designed  to  be 
substantially  larger. 

This  unusual  rotor  blade  setup  causes  various  hydrodynamic 
problems.  First  of  all,  since  a  typical  flow  incidence 
angle  near  the  hub  should  be  surprisingly  high  (e.g.,  30*), 
even  a  slight  error  in  design  may  cause  flow  separation, 
possibly  cavitation  and  then  noise  generation.  Secondly, 
even  if  design  is  made  properly,  the  same  vulnerable 
situation  is  generated  with  a  slight  flow  disturbance  or 
blade  deformation  due  to  fabrication  inaccuracy.*  A  recent 
study  at  Tetra  Tech  (see  the  report  by  Furuya,  et  al. 
(1984))  indicated  that  some  blade  deformation,  particularly 
near  the  hub,  could  cause  an  increase  of  the  power  coef¬ 
ficient,  Cp,  by  as  much  as  7.3X  (see  Figure  1-5). 
Furthermore,  there  exists  a  profound  discrepancy  between 
water  tunnel  test  results  and  actual  sea  runs.  What  causes 
such  a  discrepancy  has  not  been  clarified  to  date.  It  is 
conceivable  that  1)  a  small  trim  angle  (such  as  1  -  2*) 
existing  at  actual  sea  runs  might  have  caused  a  change  in 
boundary  layer  velocity  profile,  or  2)  the  boundary  layer 
may  be  different  between  the  water  tunnel  and  unbounded  flow 
environment  so  that  the  pumpjet  performance  is  substantially 
affected.  It  should  be  noted  that  the  utilization  of  the 
boundary  layer  is  an  advantage  in  obtaining  the  pumpjet’s 
high  efficiency  on  one  hand  but  it  is  a  disadvantage  in 
causing  many  difficult  problems  on  the  other  hand. 


*  Some  pumpjet  rotors  are  produced  by  investment  casting 
process  so  that  the  fabrication  accuracy  cannot  be 

expected  to  be  high. 
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The  turning  capability  of  the  underwater  vehicle  thrusted 
with  a  pumpjet  is  said  to  be  inferior  to  that  with,  e.g.,  a 
counter-rotating  propeller.  The  reason  for  such  inferior 
turning  capability  seems  also  attributable  to  the  utiliza¬ 
tion  of  the  wake  flow;  when  the  vehicle  turns,  the  boundary 
layer  substantially  changes.  The  pumpjet  seems  to  lose  a 
considerable  thrust  capability  due  to  the  change  of  boundary 
layer  velocity  profile,  resulting  in  a  poor  turning  capabil¬ 
ity. 

Another  problem  area  in  the  pumpjet  lies  in  the  pumpjet 
design  method.  The  only  design  method  developed  to  date  and 
used  is  a  two-dimensional  graphic  method  combined  with 
experimental  data  of  Bruce,  et  al.  (1977)  despite  the  fact 
that  the  pumpjet  experiences  a  three-dimensional  flow. 
Based  on  the  momentum  theorem  applied  to  the  cascade  con¬ 
figuration,  the  blade  sectional  pressure  increase  Ap  is 
given 


Ap  =  K  Vm'AVe  (1.1) 

where  Ap  =  local  pressure  increase  through 

the  rotor, 

Vm  =  meridional  velocity, 

AV0  =  circumferential  velocity,  and 

K  =  a  constant  determined  by  the 
cascade  configuration. 

This  two-dimensional  momentum  theory  indicates  that,  in 
order  to  generate  a  certain  pressure  increase  at  a  blade 
section,  only  the  amount  of  total  flow  deflection  in  the 
circumferential  direction  (between  the  inlet  and  exit) 
counts,  see  Figure  1-6.  In  this  method  once  the  sectional 
blade  leading  edge  and  trailing  edge  angles  are  determined, 
then  the  rest  of  the  blade  section  can  be  arbitrarily  deter¬ 
mined  by  connecting  these  predetermined  leading  and  trailing 
edges,  e.g.,  anglewise  smoothly. 


One  of  the  obvious  problems  in  this  graphic  method  arises 
from  the  fact  that  the  flow  coming  into  the  cascade  blade 
cannot  exactly  follow  the  blade  camber,  but  substantially 
deviates  from  it.  What  is  required  therefore  is  a  camber 
correction,  the  amount  of  which  depends  upon  the  cascade 
geometry.  Unfortunately,  a  typical  pumpjet  solidity*  near 
the  hub  is  larger  than  2.0  and  therefore  the  camber  correc¬ 
tion  required  there  becomes  as  much  as  5  times  in  terms  of 
lift  coefficient.  It  means  that  the  camber  graphically 
constructed  should  be  deformed  until  the  lift  coefficient 
increases  by  5  times  that  graphically  obtained.  This 
correction  is  made  semi-empirically  based  on  limited  numbers 
of  existing  experimental  data  for  cascade  blade  rlows.  In 
this  sense,  therefore,  this  graphical  method  is  useless  for 
the  blade  design  near  the  hub  and  it  can  be  said  that  the 
final  design  is  almost  entirely  dependent  upon  these  empiri¬ 
cal  data. 
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OBJECTIVES 


The  objectives  of  the  work  to  be  conducted  under  the  GHR 
program  are  therefore: 

1)  to  develop  a  more  reliable  and  accurate  pumpjet 
design  method  based  on  a  three-dimensional  pump  or 
propeller  design  theory  and  then 

2)  to  improve  the  pumpjet  performance  characteristics. 
The  characteristics  to  be  improved  include: 

a)  the  susceptibility  to  flow  disturbance  and  rotor’s 
deformation  due  to  fabrication  inaccuracy, 

b)  the  discrepancy  problem  between  the  water  tunnel 
test  results  and  high  speed  sea  runs  and 


c)  the  poor  turning  capability. 


3.0 


METHODOLOGY  SELECTION  FOR  A  THREE-DIMENSIONAL 
PUMP JET  DESIGN 


3.1  REQUIREMENTS  AND  CANDIDATES 

There  exist  several  possible  approaches  which  can  incor¬ 
porate  three-dimensionality  into  pumpjet  design  procedure. 
However,  the  following  aspects  should  be  considered  in 
selecting  such  a  methodology: 

1)  Moderate  three-dimensionality 

A  pumpjet  is  usually  installed  astern  of  under¬ 
water  vehicle  hull  where  the  hull  shape  has  a 
negative  slope  of  tapering  shape.  Although  this 
provides  three-dimensional  flow  characteristics, 
its  three-dimensionality  is  rather  mild,  unlike 
that  in  radial  pump  cases. 

2)  Capability  of  determining  detailed  blade  profile 
shapes  as  well  as  pressure  distribution 

In  the  previous  two-dimensional  graphic  method 
Bruce,  et  al.  (1974),  the  blade  profile  shape  was 
graphically  determined  for  meeting  the  head 
generation  requirement.  It  is  mainly  for  this 
reason  that  the  method  failed  to  check  the  possi¬ 
bility  of  flow  separation  after  the  blade  was 
designed.  A  new  method  to  be  developed  in  this 
research  work  should  be  the  one  with  which  the 
detailed  pressure  distribution  or  velocity  distri¬ 
bution  on  the  blade  can  be  determined. 

3)  Accurate  loading  determination  supported  by 
experiments 

When  the  sectional  loading  is  determined  analyti¬ 
cally  in  the  course  of  designing  a  pumpjet,  it  is 
usually  quite  inaccurate  since  such  loading 
substantially  changes  due  to  the  effect  of  adja¬ 
cent  blades.  It  is  therefore  necessary  for  the 
new  method  to  incorporate  the  cascade  effect  into 
design  procedure,  or  to  use  an  empirical  approach 
to  increase  such  accuracy. 

With  these  features  taken  into  consideration,  the  following 
three  candidate  methods  are  compared  in  Section  3.2: 

Method  I:  Katsanis*  Quasi-Three-Dimensional  Method 

Method  II:  Blade-Through  Flow  with  Blade-to-Blade 

Flow  Method 


Method  III:  Singularity  Distribution  Method 


for  which  simple  explanation  will  be  given  in  the  following. 

Method  I:  Katsanis*  (1964)  Quasi-Three-Dimensional  Method 

In  this  method  it  is  first  assumed  that  a  mean  stream 
surface  from  hub  to  shroud  between  blades  is  known  in 
advance.  On  this  stream  surface  a  two-dimensional  solu¬ 
tion  for  the  velocity  and  pressure  distributions  is 
obtained.  Then,  an  approximate  calculation  of  the  blade 
surface  velocities  is  made.  This  method  is  based  on  an 
equation  for  the  velocity  gradient  along  an  arbitrary 
quasi-orthogonal  rather  than  the  normal  to  the 
streamline.  Since  the  solution  is  obtained  on  this 
quasi-orthogonal  line,  in  this  method*  an  iteration  pro¬ 
cedure  needed  in  the  previous  orthogonal-line  methods 
can  be  eliminated  and  a  solution  can  be  obtained  in  a 
single  computer  run. 


Method  II:  Blade-Through  Flow 
Method 

with 

Blade-to-Blade 

Flow 

The  blade-through  flow  is 

first 

obtained  by, 

e .  g  • , 

Streamline  Curvature  Method 

( SCM ) . 

Once  the  stream 

sur- 

face  is  found,  it  is  mapped 

to  a 

two-dimensional 

plane 

so  that  the  blades  cut  through  by  the  stream  surface 
become  a  row  of  blades,  i.e.,  cascade  on  a  plane.  On 
this  cascade  configuration,  the  blade-to-blade  flow  will 
be  solved.  Difficulty  in  doing  this  lies  in  the  fact 
that  the  governing  equation  is  not  a  Laplace  equation 
any  more  on  this  two-dimensional  plane,  but  a  Poisson 
equation  due  to  the  deviation  of  stream  surface  from  a 
perfect  cylinder.  In  order  to  account  for  such 
deviation  into  the  two-dimensional  flow,  appropriate 
source/sink  and  vortices  should  be  distributed  over  the 
entire  flow  field.  This,  in  turn,  results  in  the  change 
of  blade  camber  shape.  The  design  procedure  relies  on 
an  iteration  scheme. 


Method  III*  Singularity  Distribution  Method 


The  method  is  similar  to  that  used  in  design  of  conven¬ 
tional  propellers,  see,  e.g.,  the  work  by  Kerwin  and 
Leopold  (1964).  The  differences  in  velocities  between 
the  pressure  and  suction  sides  of  rotor  blade  can  be 
represented  by  distributions  of  singularities  such  as 
source/sink  and  vortex.  The  strengths  of  such  singu¬ 
larities  are  determined  by  satisfying  the  boundary 
conditions  on  the  blade  surface  as  well  as  those  at 
infinity.  The  methodology  is  described  in  the  paper  of 
Kerwin  and  Leopold  (1964)  in  detail. 

The  disadvantage  of  the  method  lies  in  the  computational 
complexity  and  instability.  Furthermore,  this  type  of 
method  is  suitable  for  design  of  devices  used  in  the 
open  field,  but  not  so  for  those  used  in  the  internal 
flow  since  it  does  not  take  advantage  of  confined  flow 
configuration  available  for  the  latter  case. 

3.2  COMPARISON  AND  SELECTION 

Table  3-1  provides  qualitative  evaluation  on  three  candidate 
methods  described  in  the  previous  section  over  various 
hydrodynamic,  numerical  and  design  aspects.  As  seen  from 
this  table,  a  combination  of  blade-through  method  with 
blade-to-blade  flow  seems  to  have  an  advantage  over  the 
other  two  methods.  Particularly,  the  method  has  the  capa¬ 
bility  of  determining  detailed  blade  profile  shape  as  well 
as  loading  and  velocity/pressure  distribution  with  accuracy 
verified  by  existing  cascade  experimental  data.  It  is  for 
this  reason  that  the  blade-through  flow  with  blade-to-blade 
flow  method  has  been  selected  as  a  basic  concept  for  devel¬ 
oping  a  three-dimensional  pumpjet  design  method. 

Similar  methods  already  exist  for  design  of  quasi-axial 
pumps  and  compressors.  However,  those  methods  have  many 
inadequate  features  in  their  design  procedure.  Furthermore, 


it  is  assumed  in  these  design  methods  that  the  incoming  flow 
is  more  or  less  uniform,  unlike  the  pumpjet  where  the  rotor 
should  be  designed  for  highly  retarded  velocity  distribution 
due  to  viscous  boundary  layer  on  the  hull.  The  following 
section  describes  the  blade-through  flow  (BT)  with  blade-to- 
blade  flow  (BTB)  method  with  various  aspects  of  modifica¬ 
tions  and  improvement  necessary  for  developing  the  pumpjet 
design  method. 


4.0 


SELECTED  DESIGN  METHOD— BLADE-THROUGH 
BLADE-TO-BLADE  METHOD 


Design  of  a  pumpjet  for  an  underwater  vehicle  requires  pre¬ 
liminary  information  on  the  vehicle  including  its  geometry 
and  hydrodynamic  drag  coefficient.  Furthermore,  most 
importantly,  the  velocity  profile  at  an  upstream  reference 
section  should  be  obtained  either  analytically  or 
experimentally.  Any  error  in  the  velocity  profile  would 
result  in  a  pumpjet  of  lower  efficiency  or  failure  of  the 
pumpjet  meeting  the  specifications  at  the  design  point.  In 
the  present  study,  it  is  assumed  that  this  velocity  profile 
is  given  at  a  goal  speed  or  at  the  corresponding  Reynolds 
number . 

The  first  step  for  design  of  a  pumpjet  (see  Figure  4-1)  is 
to  determine  the  shroud  intake  diameter.  From  the  viewpoint 
of  cavitation,  the  maximum  and  minimum  shroud  diameter  to 
prevent  cavitation  must  exist.  If  it  is  too  large,  the 
rotor  blade  tip  speed  becomes  too  high  so  that  cavitation 
occurs.  On  the  other  hand,  if  it  is  too  small,  the  rotation 
speed  must  be  increased  to  generate  the  required  head  so 
that  the  chance  of  cavitation  inception  also  increases. 
Another  aspect  of  determining  the  shroud  diameter  stems  from 
consideration  of  the  overall  propulsive  efficiency.  The 
equation  for  global  momentum  balance  should  be  able  to 
determine  an  efficiency-optimum  shroud  diameter  for  the 
given  velocity  profile  and  vehicle  drag.  The  detailed 
mathematical  formulation  and  sample  calculations  will  be 
given  in  Section  4.1. 

Once  the  shroud  diameter  is  determined,  the  streamline  will 
be  calculated  by  using  the  streamline  curvature  method 
(SCM).  In  this  calculation,  the  loading  distribution  on  the 
rotor  and  blade  thickness  must  be  assumed  in  advance.  One 
of  the  major  concerns  in  using  the  existing  streamline  cur¬ 
vature  method  lies  in  the  fact  that  SCM  may  only  be  used  for 
relatively  uniform  incoming  flow,  but  may  generate  a  sub- 


stantial  error  for  a  thick  wake  flow,  i.e.,  highly  retarded 
flow  due  to  the  viscous  boundary  layer  on  the  vehicle  hull. 
Detailed  mathematical  formulation  and  sample  calculations 
are  presented  in  Section  4.2.  Also  included  are  discussions 
regarding  the  problems  of  application  of  conventional  SCM  to 
the  thick  wake  flow. 

The  next  step  of  the  design  method  is  to  map  the  stream  tube 
or  surface  calculated  by  SCM  onto  a  plane  so  that  the  rotor 
blades  are  mapped  into  cascade  configuration.  If  the  stream 
surface  is  totally  cylindrical  shape,  the  governing  equation 
to  be  used  for  the  cascade  analysis  will  be  a  Laplace 
equation.  Unfortunately,  the  stream  surface  is  of  three- 
dimensional  cone  shape  in  general  for  the  tail  cone  section 
of  the  underwater  vehicle.  The  field  governing  equation  now 
becomes  a  Poisson  equation,  for  which  the  results  of  power¬ 
ful  potential  theory  analysis  are  no  more  applicable.  A 
method  of  correcting  the  effect  of  the  Poisson  equation  on 
the  potential  theory  results  is  introduced  to  modify  the 
blade  profile  shape  obtained  in  the  potential  theory.  In 
choosing  the  blade  profile  shape,  the  experimental  data  are 
used  to  ensure  that  there  is  no  chance  of  flow  separation 
due  to  overloading  on  the  blade.  Furthermore,  based  on  the 
calculated  velocity  along  the  blade,  the  cavitation  incep¬ 
tion  is  checked.  If  there  exists  a  chance  of  either  flow 
separation  or  cavitation,  the  loading  distribution  from  hub 
to  tip  should  be  changed.  If  such  a  change  is  made,  and/or 
thickness  of  blades  is  changed,  the  streamline  curvature 
method  should  be  used  again  to  determine  the  new  location  of 
streamline  or  stream  surface.  This  iterative  procedure  will 
be  repeated  until  an  overall  convergent  solution  is 
obtained.  Section  4.3  describes  the  technical  approach  to 
be  used  for  the  blade-to-blade  flow  analysis. 
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PUMP JET  GLOBAL  HYDRODYNAMICS 


It  is  a  well-known  fact  that  the  pumpjet  utilizes  the  tail¬ 
cone  low-energy,  boundary  layer  flow  in  order  to  achieve  its 
high  efficiency.  It  means  that  the  optimum*  pumpjet  design 
depends  entirely  upon  the  incoming  flow  velocity  profile. 


Figure  4-2  shows  a  schematic  diagram  of  an  underwater 
vehicle  tail  cone/pumpjet  flow.  (7)  and  (7)  in  Figure  4-2 
are  considered  to  be  the  upstream  and  downstream  reference 
stations,  respectively,  where  it  is  assumed  that  the  free- 
stream  static  pressure  exists,  whereas  (7)  and  (T)  are  the 
rotor  inlet  and  exit  stations. 


4.1.1 


Calculation  of  Thrust  Force 


The  thrust  force,  T,  due  to  the  pumpjet  work  can  be  deter¬ 
mined  by  applying  the  momentum  equation  to  the  control 
volume  enclosed  by  stations  (7) ,  (7)  and  the  stagnation 
streamline  (see  Figure  4-2). 


T  =  J  dm7« V7(r  )cose7  -  J  dm^V^r  Jcossj^ 

r  1 1  -»  fiii 


(4. 1. 1-1) 


where  dm7  =  p2irrV  7  (r)  cose  7 


(4. 1. 1-2) 


1^  =  p2TrrV  ^  ( r )  cose  ^ 


V.(r),  V7(r)  =  meridional  flow  velocities  at  (T)  and 

(7),  respectively. 

e.,  e7  =  meridional  flow  angles  at  (7)  and 
(7),  respectively. 


Note 


*By  w optimum**  pumpjet  design  it  means  that  of  proving 
the  maximum  propulsion  efficiency. 
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Defining  the  following  quantities, 
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and  normalizing  T,  Eqn.  (4. 1.1-3)  becomes 
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(4. 1.1-6) 


where  the  mass  conservation  equation  below  has  been  used; 
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(4. 1. 1-7) 
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then  we  obtain,  from  (4. 1.1-6) 
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4.1.2  Relationship  Between  Pump  Head  and  Thrust 
The  hydraulic  head  of  the  pump  Hr  is  given 
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where  VjCr),  VgjCr)  =  meridional  and  circumferential  veloci¬ 
ties  at  station  (3),  respectively. 


With  the  hydraulic  efficiency  r»R  introduced,  the  actual  head 
generated  in  the  fluid  is  H, 

H  =  nR*HR  (4. 1.2-2) 

where  H  can  be  defined 
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Ki  =  h^ad  loss  coefficient  between  station  (T)  and 
(?)  (but  mostly  inlet  loss  and  see  Appendix 
for  determining  K^). 


Since 


AV  =  V  7  -  V, 
m  7  1 


(4. 1.2-6) 


and  AV m  can  be  approximated  as  AV m  -  V7  -  V  ^  Eqn. 


(4. 1.2-3)  becomes 
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From  Eqns .  (4.1.1-12)  and  (4. 1.2-7), 
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4.1.3  Power  Calculation 


The  power  to  be  used  on  the  pumpjet  rotor  shaft  can  be 
calculated  by  integrating  the  product  of  the  local  head  and 
the  local  mass  flow  rate  over  the  entire  duct  flow, 
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From  (4. 1.2-1), 
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Therefore,  the  propulsive  efficiency  np  can  be  calculated 
from  the  following  definition; 

CT 


nP  =  Cr 


(4. 1.3-3) 


4.1.4  Procedure  of  Calculating  Propulsive  Efficiency  with 
Vi(r),  etc,  given  and  Sample  Results 


Now  all  tools  for  calculating  the  propulsive  efficiency  rip 
with  Vi(r)  given  are  provided.  The  principal  equations  to 
be  used  will  be  Eqns.  (4.1.1-13),  (4. 1.2-6),  (4. 1.2-8), 


(4. 1.3-2),  and  (4. 1.3-3).  A  flow  chart  describing  the 
calculation  procedure  is  given  in  Figure  4-3. 

Sample  velocity  profiles,  V^(r)  and  V2(r),  are  shown  in 
Figures  4-4  and  4-5.  The  flow  angles,  9]_,  92  and  97,  can  be 
obtained  from  the  drawings  of  a  typical  underwater  vehicle 
tail  cone  profile.  The  head  loss  coefficient,  ,  can  be 
calculated  from  the  formula  given  in  Appendix  A  with  the 
velocity  distributions  and  pressure  distributions  both  at 
stations  (T)  and  (2),  which  are  also  given  in  Figures  4-4 
and  4-5. 

Figures  4-6  to  4-7  show  the  results  of  calculations  made  for 
various  parameters,  including  the  shroud  opening  diameter, 
r^/rg,  the  exit  jet  flow  angle  97,  and  flow  velocity  ampli¬ 
fication  factor.  The  design  values  of  r^/rg  and  flow  ampli¬ 
fication  factor  are  known,  i.e., 

( r l/rg )g  =  .93. 

Flow  Amplification  Factor  =  1 

but  that  of  97  is  not  known  except  for  the  fact  that  the 
average  geometric  angle  of  the  shroud  and  tail  cone  angle  at 
exit  is  about  11*.  Theoretically,  however,  the  jet  coming 
out  from  the  shroud  exit  with  11*  should  align  itself  in  the 
direction  of  the  body  axis,  indicating  that  97  could  be 
zero.  The  present  calculations  were  therefore  made  for 
97  =  O’ ,  5* ,  8’  and  11* . 

The  shroud  opening,  r^/rg,  was  also  varied  in  the  present 
analysis  in  order  to  determine  the  optimum  shroud  opening 
radius  in  terms  of  efficiency.  Also  changed  was  the 
incoming  flow  velocity  amplitude  to  investigate  a  possibil¬ 
ity  of  pumpjet  efficiency  improvement  in  combination  with 
the  tail  cone  flow  pattern  change. 

The  pump  hydraulic  efficiency  tir  for  this  type  of  pump,  used 
in  the  analysis,  is  about  89*,  which  is  the  measured  value 
by  many  pump  makers. 


Figures  4-6  to  4-8  show  the  calculated  propulsive  efficien¬ 
cies  as  a  function  of  the  shroud  opening  r^/rg  with  the  flow 
amplification  factor  as  a  parameter  for  07  =  0*,  5*  and  11*, 
respectively.  As  can  be  seen  from  these  figures,  the  effi¬ 
ciency  curve  has  the  maximum  value  at  an  r^/rg  value  spe¬ 
cific  for  the  conditions  used. 

Figure  4-6  shows  that  the  efficiency  at  the  design  condition 
should  be  88.4%  when  07  =  0*  is  assumed.  The  design  shroud 
opening,  ri/rg  =  .93,  is  slightly  on  the  smaller  side  than 
that  for  the  maximum  efficiency.  The  maximum  efficiency  of 
88.8%  can  be  obtained  at  a  slightly  larger  shroud  radius, 
i.e.,  r^/rg  =  .95.  Also  seen  from  Figure  4-6  is  the  fact 
that  the  smaller  the  incoming  flow  velocity,  the  larger  the 
maximum  propulsive  efficiency.  This  indicates  that,  if  the 
incoming  velocity  amplitude  at  the  tail  cone  area  can  be 
reduced  by  some  means,  the  propulsive  efficiency  of  the 
pumpjet  is  substantially  increased. 

It  should  also  be  remembered  that  the  actual  efficiency 
achieved  is  76.9%,  much  lower  than  any  of  the  values  calcu¬ 
lated  here. 

As  07  increases,  the  efficiency  curves  shift  to  the  lower 
efficiency  side,  see  Figures  4-6  to  4-8.  This  is  naturally 
expected  since  the  jet  thrusting  force  is  not  effectively 
utilized  as  07  increases. 

It  may  also  be  coincidental,  as  seen  in  Figure  4-8,  that  if 
07  =  11*  is  used  as  obtained  from  the  pumpjet  exit  geometry, 
the  current  shroud  opening,  r^/rg  =  .93,  is  the  optimum 
selection  for  providing  the  maximum  efficiency,  82.9%, 
smaller  than  that  of  the  original  design.  It  should  be 
pointed  out  that  if  07  =  11*  is  the  true  exit  jet  flow 
angle,  an  increase  in  the  shroud  radius  for  alleviating  the 
flow  separation  problem  may  cost  a  substantial  efficiency 
loss  (see  Figure  4-8).  On  the  other  hand,  if  07  =  0*  is  the 
true  value,  a  moderate  increase  (e.g.,  2%)  in  the  shroud 


radius  will  increase  the  efficiency  in  addition  to  the  effi¬ 
ciency  gained  due  to  the  suppression  of  the  flow  separation. 
However,  again  the  penalty  exists  when  the  amount  of  shroud 
radius  increase  exceeds  more  than  4%. 

Figure  4-9  summarizes  the  present  calculations  in  terms  of 
the  maximum  propulsive  efficiency  with  07  and  the  flow 
amplification  factor  as  parameters.  It  is  shown  that  a 
substantial  efficiency  improvement  may  be  achieved  by 

1)  reducing  the  incoming  flow  velocity  amplitude 
by  modifying  the  tail  cone  profile  shape  and 
thus  changing  the  boundary  layer  flow, 

2)  choosing  the  optimum  shroud  opening  radius 
depending  on  the  flow  conditions  (see  also 
Figures  4-6  to  4-8), 

3)  redirecting  the  jet  flow  at  the  shroud  exit  as 
close  to  the  body  axis  as  possible,  if  the 
jet  flow  of  the  current  design  is  not. 

4.2  BLADE-THROUGH  FLOW  ANALYSIS  -  STREAMLINE  CURVATURE 

METHOD  ( SCM ) 

4.2.1  Mathematical  Formulation 

From  the  definition  of  entropy,  S,  in  the  second  law  of 
thermodynamics,  the  following  relationship  is  obtained  for  a 
reversible  transformation,  i.e., 

TVS  =  VQ  (4.2. 1-1) 

where  T  is  the  temperature  and  VQ  is  the  amount  of  heat  the 
system  under  consideration  receives.  On  the  other  hand,  the 
first  law  of  thermodynamics  says 

VE  =  VQ  +  VW  (4.2. 1-2) 

where  E  is  the  internal  energy  and  W  is  the  work  performed 
on  the  fluid.  Since  VW  =  -pVv,  (4.2. 1-2)  becomes 


=  VQ  -  p7(±) 


(4.2. 1-3) 


i 
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B 


[ 
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where  v  is  the  specific  volume  of  fluid  and,  in  terms  of 
fluid  density,  p,  v  =  1/p.  The  definition  of  enthalpy,  H, 
is  given  by 


H=^u2  +  E+  £  +  |)  (4.2. 1-4) 

where  u  is  the  amplitude  of  flow  velocity  and  ip  is  the 
potential  energy.  Gradient  of  H  yields 


7H  =  7(^  u2  +  Ji)  +  7E  +  Vp  +  p7(£) 

From  the  above  equation  and  Eqn.  (4. 2. 1-3), 

7H  =  T7S  +  7(^  u2  +  V)  +  ±  Vp  (4.2. 1-5) 

The  steady-state  momentum  theorem  gives 

pu  •  7u  =  -7p  -  p  7$J  (4.2. 1-6) 

where  an  assumption  has  been  made  that  -7<J  *  F,  where  F  is 

an  external  force.  By  using  a  vector  identity,  _u  x  (7  x  u) 
-  1  7  u2  -  tj  •  7lj,  Eqn.  (4.2. 1-6)  is  now  written 

1  2  i 

ux  (7xu)  =  ^  7u  +  7{J  +  -i  7p 

2 

or  u  x  w  =  7  (y-  +  V)  +  “  7p  (4.2. 1-7) 

where  u  =  7  x  u.  Substituting  Eqn.  (4.2. 1-5)  into  (4.2. 1-7) 
gives 


u  x  0)  =  7H  -  T7S ,  (4.2. 1-8) 

a  relation  first  found  by  Crocco  (1937),  which  will  be  used 
to  derive  the  formula  used  for  SCM  hereafter. 

By  using  the  cylindrical  coordinate  system  (r,e,z),  the 
velocity  components  of  u  are  defined  by 

u  =  (ur,  Uq ,  uz) .  (4.2. 1.9) 
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Thus,  the  components  of  vortex  term  w  is  written 


(a)  wr 
Ch)  o>e 


(V  x  u)r 
(V  x  u)0 


(4.2.1-10) 


By  introducing  a 
•  4-10) 


=  (V  x  u)z 
direction 


-  I  /*rU9 
r  \  3r 

m,  defined 


9uA 
3e  J 

by,  (see  also  Fig. 
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dm  t  dr 
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um  :  u„  :  u, 
m  r  z 

1 

2  2 
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• 

(b) 

U  =  U 

m  r 

+  U2 

(c) 

tan  = 

ur/uz 

(4.2.1-11) 

(d) 

u  =  u 
r  m 

simj) 

• 

(e) 

u  =  u 
z  m 

COS(|)  , 

it  becomes  evident  that  the 

m-direction  is  the  "meridional 

to 

direction'*  or 

on  the 

projection  of  streamline  in  the  r-z 

plane.  The  directional  derivative  with  respect  to  m  then 

becomes 

• 

• 

(a) 

u  |—  = 
m  am 

u  I3 
r  3r 

+  u  f- 
z  3z 

(b) 

3 _ 3r 

i-  + 

3z  i_ 

I 

3m  ~  3m 

3r 

3m  3z 

=  sin$  +  cos<f>  (4.2.1-12) 

where  (4.2.1-11)  has  been  used. 

Furthermore,  the  e-component  of  the  Crocco  equation 
(4.2. 1-8)  gives 


(4.2.1-13) 
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_  J.  / 3H  as\ 
r  y3e  "  3e/ . 

Under  the  assumption  of  axisymmetricity , 


3H 

39 


=  0 
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3S 

39 


=  0 


9 


so  that  Eqn.  (4.2. 1.-13)  becomes 

3u  3u  3ru  3ru 

ur  3©“  +  uz  3e  =  ^  Tr  +  ^  ~3z 


or 
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3_ 

39 


(u2  + 


3rU8 
um  3m 


Using  (4.2.1-llb), 

3  u  3rua 

_ m  _9 

39  ■  3m 


(4.2.1-14) 


It  is  now  ready  to  perform  a  coordinate  transform  of  Eqn. 
(4.2.1-10)  by  using  Eqns.  (4.2.1-11),  (4.2.1-12)  and 
(4.2.1-14).  The  result  is 


(a)  (£r  =  (Vxu)  r 


(b)  wQ  =  (Vxu)0 


.  l  f3ru  „  3ru  _ 

j-j-E  .  sin* 

cos*  ft] 

1  [  2  /simt)  3um 

[  m  \  um  3m 

(4.2.1-15) 

cosfA  1  3um] 

rm  /  ”  2  3r  J 


(c)  wz  =  (Vxu)z 
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3m 


where  rm  is  the  radius  of  curvature  of  the  meridional 
streamline  projection,  defined  by 
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(4.2.1-16) 


1 _ 3^ 

rm  "  "  3m 
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Further  application  of  axisymmetr icity  to  Eqn.  (4.2.1-15) 
yields 
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where  the  following  relations 


It  .  o 

3e  "  0 


3ru„  3u 

_ 9  _  _ m 

3m  ~  3e 
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(4.2.1-17) 


have  been  applied. 


The  right-hand  side  of  Eqn.  (4.2. 1-8)  becomes  VH  under  the 
assumption  of  adiabatic  process  for  the  fluid  to  go  through 
the  pump  channel. 

Now,  it  is  ready  to  write  the  r-component  of  Eqn.  (4.2. 1-8) 
for  the  meridional  flow  velocity  and  the  final  form  is  shown 
after  some  rearrangement: 


2 


(■ 


sinij)  3u  m 

um  3m 
m 


(4.2.1-18) 


Eqn.  (4.2.1-18)  can  he  written  as 


6 


25 


6 


» 

i 

II 


P(  r ) 


m 
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where 


P( r )  =  2 


TCr)  =  2 


3H  u.  3rue\ 
3r  r  3r  / 


(4.2.1-19) 


(4.2. 1-20 ) 


In  Eqns.  (4.2.1-18)  -  (4.2.1-20),  the  first  term  of  P(r), 
i.e.,  sin^)/um  •  3um/3m,  will  provide  some  difficulty  in 

numerical  computations  since  it  is  related  to  the  deriva¬ 
tives  with  respect  to  *'mH.  The  basic  philosophy  of  the 
streamline  curvature  method  is  to  express  the  meridional 
velocity  in  terms  of  “r"  and  "r-derivatives"  so  that  um  can 
be  solved  in  the  direction  of  r  only.  This  feature  will  be 
of  advantage  in  numerical  computations  since  the  derivatives 
with  respect  to  MmH  are  not  needed  and  thus  the  m- 
directional  control  points  do  not  have  to  be  taken  in  fine 
increments.  Fortunately,  sin^/um  •  3um/3m  can  be  expressed 
in  terms  of  r  by  using  the  continuity  equation 


1 

r  3r 


+ 1 

r  3e 


0. 


Again,  applying  the  axisymmetric  assumption  and  Eqns. 
(4.2.1-11),  (4.2.1-12)  and  (4.2.1-16),  the  following  rela¬ 
tion  is  obtained 


sin<fr  3um 
um  '  9m 


sin* 

r 


(1  + 


m 


_r _ 

cos  j) 


)  +  tan  $>•  || 


(4.2.1-21) 

The  basic  equation  for  the  streamline  curvature  method, 
i.e.,  Eqn .  (4.2.1-18),  is  expressed  in  terms  of  wr"  except 
for  the  radius  of  curvature,  rm,  so  that  it  can  be  readily 
solved  numerically.  The  only  problem  remaining  is  that 
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um(r)  cannot  be  uniquely  determined.  This  problem  can  be 
resolved  by  applying  the  mass  conservation  equation 


(4.2.1-22) 


where  K&  is  the  blockage  factor  due  to  the  blade  displace¬ 
ment  thickness  as  well  as  that  of  the  boundary  layer,  and  ^)q 
is  the  angle  between  the  line  of  integration  (called 
“q-line"  hereafter)  and  the  line  normal  to  the  streamline. 
With  this,  the  mathematical  formulation  for  the  streamline 
curvature  method  (SCM)  is  completed.  In  what  follows,  the 
numerical  solution  method  for  SCM  will  be  described  in 
detail . 


4.2.2  Solution  Method 

Since  Eqns.  (4.2.1-18)  and  (4.2.1-22)  are  highly  nonlinear 
for  um,  only  an  iterative  procedure  depending  on  numerical 
analysis  is  a  possible  solution  method.  First  of  all,  it  is 
assumed  that  the  distribution  of  upstream  flow  velocity  is 
known  as  function  of  r.  Figure  4-11  shows  a  sample  flow 
configuration  on  an  underwater  vehicle  tail  cone  where  the 
upstream  location  in  this  case  is  identified  by  station  1  (I 
=  1).  The  upstream  flow  velocity  can  either  be  uniform  or 
nonuniform*.  Station  1  is  then  divided  into  a  finite  number 
of  control  points  including  the  hub  and  the  inside  wall  of 
the  shroud.  In  Figure  4-11,  a  total  of  6  control  points 
(J  =  6)  are  used.  By  using  the  mass  conservation  equation 
(4.2.1-22),  the  local  mass  flow  rate  (gj,j+i)  between  each 
two  adjacent  control  points  J  and  J+l  is  calculated.  The 
total  mass  flow  rate,  Gj,  is  just  the  summation  of  local 
flow  rates 


*Note )  However,  the  upstream  flow,  which  is  severely 
retarded  or  highly  nonuniform  due  to,  e.g.,  the 
viscous  effect,  may  present  problems  of  accuracy, 
which  will  be  discussed  in  Section  4.2.3. 


=  j,  j+i 


(4. 2.2-1) 
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G 


J+l 


where  9lfJ,J+l  2lr  /  Kb‘pr  um1(r)  cos^qdr  (4. 2. 2-2) 


For  Stations  I  =  2,3,4...,  the  initial  control  points  and 
initial  um  velocity  profile  are  also  needed  for  the  itera¬ 
tion  procedure.  For  the  case  of  handling  a  uniform  upstream 
flow  velocity,  the  selection  of  control  points,  J  =  1  ~  6, 
and  determination  of  initial  flow  velocity  may  be  done  in 
the  same  manner  as  that  for  the  first  Station,  1=1,  since 
the  constant  velocity  distribution  can  be  assumed.  However, 
the  case  of  nonuniform  flow  velocity  distribution  will 
require  a  little  care  for  selection  of  control  points  and 
determination  of  initial  flow  velocity  distribution.  Among 
many  possible  ways,  it  has  been  decided  herein  that  the 
velocity  distribution  is  assumed  to  have  a  similarity  nature 
as  that  of  the  upstream  at  I  =  1 ,  i.e., 


uml(r)  = 


kIu«Utr) 
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r-r 


r  HI  +  r 


HI 


Sl“rHI 


(rSl"rHl) 


I  =  2,3. . . 


where 


(4. 2. 2-3) 


r_,  r_.  =  radius  of  shroud  internal  wall  at  Station  1 
’  and  I  (>  2),  respectively 

rHl  rHI  =  radius  hub  at  Station  1  and  I  (>  2), 

’  respectively 

kj  =  an  arbitrary  constant,  dependent  of  Station  I, 
to  be  determined  later. 


An  arbitrary  constant  kj  is  used  to  adjust  the  total  flow 

• 

rate  at  Station  I  >  2  becomes  Gi  when  umi(r)  is  substituted 
into  Eqn.  (4.2.1-22).  Once  kj  is  properly  determined,  the 
control  points  J  =  2,3...  can  be  determined  by  using  Eqn. 
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(4. 2. 2-2).  In  order  to  carry  out  the  above  computations, 
should  be  known  in  advance.  If  the  control  points  at  every 
control  station  are  known,  ()  can  be  calculated  by  connecting 
these  points  for  each  streamline  by  using,  e.g.,  cubic 
spline  method.  However,  at  the  first  iteration,  even  these 
points  are  not  yet  known.  Therefore,  £  should  be  determined 
by  guess.  One  possible  way  is  to  linearly  interpolate  ^  for 
2  <  J  <  5  from  (>  at  0  =  1  and  £  at  J  =  6,  i.e.,  the  hub  wall 
angle  and  shroud  internal  wall  angle,  respectively. 

It  is  now  ready  to  calculate  a  new  set  of  um’s  or  3um/3r  at 
I  =  2,3,4...  by  using  Eqn.  (4.2.1-18).  Since  actual  calcu¬ 
lations  are  made  on  3um/3r,  an  integration  constant  should 
be  determined  to  uniquely  determine  um  itself.  This 
constant  can  be  readily  determined  by  applying  the  mass  con¬ 
servation  equation  (4.2.1-22),  the  control  points  ( 3  = 

2,3...)  at  each  I  >  2  should  be  shifted  according  to  Eqn. 
(4. 2. 2-2).  Integral  limits,  rj’s,  which  determine  the 
control  points  on  the  q-line  for  3  =  2,3...  are  determined 
one  by  one  starting  from  the  hub  in  such  a  way  that  the  mass 
flow  rate  gi,j,j+i  remains  the  same  in  each  stream  sheet  as 
that  for  1  =  1.  This  iteration  process  must  be  repeated 
until  convergence  for  um’s  as  well  as  the  location  of 
control  points  (or  streamlines)  is  obtained. 


4.2.3  Potential  Problems  of  SCM  for  Highly  Nonuniform 
Velocity  Profile  Due  to  Viscosity 

The  present  streamline  curvature  method  (SCM)  to  be  used  for 
determining  the  meridional  flow  streamlines  is  based  on  the 
momentum  equation  and  mass  conservation  equation  with 
viscous  effects  totally  ignored.  Therefore,  if  the  upstream 
flow  is  the  one  fully  retarded  due  to  the  viscous  effect, 
i.e.,  boundary  layer  flow,  so  that  the  velocity  distribution 
is  highly  nonlinear,  the  application  of  the  present  SCM  may 
create  substantial  inaccuracy  in  determining  the  location 
and  velocity  of  streamlines.  This  point  is  clearly 
understood  by  investigating  the  momentum  equation 


(4.2.1-19)}  P(r)  and  T(r)  are  only  dependent  upon  r  except 
for  rm  which  is  a  function  of  curvature  of  streamline.  It 
means  that  um  is  a  weak  function  of  the  axial-direction 
coordinate  so  that  um  at  a  certain  station  is  almost 
entirely  determined  by  the  inner  and  outer  wall  curvature. 
No  matter  how  strongly  the  incoming  flow  velocity  is 
retarded,  the  flow  velocity  will  become  more  or  less  uniform 
before  the  flow  travels  too  far  downstream  because  the  cur¬ 
vature  effect  (rm)  cannot  last  too  long. 

The  above  discussions  seem  to  suggest  that  the  development 
of  streamline  curvature  method  (SCM)  with  viscous  effect 
incorporated  may  be  in  order,  particularly  for  handling  the 
highly  viscous  flow  near  the  tail  cone  area  of  underwater 
vehicle . 

The  momentum  equation  for  such  a  flow  should  be  of  the  form 
u  x  u  =  VH  -  TVS  -  vV2u  (4. 2. 3-1) 

instead  of  Eqn.  (4. 2. 1-8). 

One  additional  term  will  make  the  problem  extremely  complex 
and  this  problem  will  be  handled  in  the  FY-86  GHR  program. 

4.2.4  Numerical  Results 

The  streamline  curvature  method  described  in  Sections  4.2.1 
and  4.2.2  was  used  to  calculate  the  streamlines  for  a  typi¬ 
cal  underwater  vehicle  tail  cone  area  with  a  shroud.  In  the 
present  case  the  upstream  flow  velocity  was  assumed  to  be 
uniform.  A  total  of  14  q-lines  (I  =  14)  were  used  with  5 
control  points  (J  =  5)  at  each  q-line,  see  Figure  4-12.  The 
solid  lines  are  the  initial  guess  for  the  streamlines 
whereas  the  dashed  lines  are  the  converged  solution  for  the 
final  streamlines.  It  is  seen  from  this  figure  that  these 
two  sets  of  lines  match  well  to  each  other  except  for  the 
area  behind  the  middle  chord  of  rotor.  It  means  that  the 
initial  guess  used  here  was  very  accurate  until  the  flow 


passes  the  rotor  and  stator.  Due  to  the  initial  accurate 
guess  for  streamlines,  computer  time  was  minimal. 

4.3  BLADE-TO-BLADE  FLOW 


4.3.1  Transformation 


Under  the  assumption  that  an  axisymmetric  stream  surface 
exists  in  a  rotating  machine,  from  the  conservation  equation 
of  circulation,  i.e,  ^  x  _w  +  2w  =  0,  the  following  relation 
is  obtained  for  the  relative  flow, 


9w 
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3(rw  Q) 


=  2uir 


3r 


(4.3. 1-1) 


3©  3m  '  3m 

where  wm  and  we  are  relative  flow  velocities  in  the  direc¬ 
tion  of  m  and  e,  see  Figure  4-13.  The  continuity  equation 
for  the  same  stream  surface  is  also  written 


3(  bpw  Q) 


3( bprw  J 


=  0 


3©  3m 

where  b  is  the  thickness  of  stream  surface. 


(4.3. 1-2) 


Then,  a  stream  function  can  be  defined  as 
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(4.3. 1-3) 


Substitution  of  we  and  wm  in  Eqn.  (4.3. 1-3)  into  Eqn. 
(4.3. 1-1)  yields 


32|P  32|J  / 1  3r  1  3b\  3^ 

r23e*  3m2  l  r  3  m  ~  5”  3  m  /  3  m 


-2bpw  sinX 


(4.3. 1-4) 


where  X  is  the  angle  of  the  line  tangent  to  the  stream  sur¬ 
face  at  the  point  of  interest  made  with  the  axis  of  rota¬ 
tion,  see  Figure  4-13. 


This  three-dimensional  axiymmetric  stream  surface  can  be 
mapped  onto  a  two-dimensional  plane,  (X,  V),  see  Figure 
4-14,  by  the  following  mapping  functions 


where  r0  is  an  arbitrary  constant  which  is  used  for  the  pur¬ 
pose  of  scaling  between  the  physical  coordinate  space  and 
mapped  plane,  (X,  Y).  The  governing  equation  (A. 3. 1-4)  can 
now  be  written  in  the  (X,  Y)  coordinate  system  by  using  Eqn. 
(4.3. 1-5) 

r  2 

V2lp  =  -2bpu  (  — )  sinX 
o 
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bp)  3X  3  X  3  Y  3  Y  j . 

Also,  the  relative  velocities  in  the  X 
given 


(4.3. 1-6) 
and  Y-directions  are 
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(4.3. 1-7) 


As  seen  from  Eqn.  4.3. 1-6,  the  governing  equation  for  the 
(X,  Y)  plane  is  now  a  Poisson  equation  instead  of  the 
Laplace  equation,  which  exists  only  for  a  flow  on  a  per¬ 
fectly  cylindrical  stream  surface.  Therefore,  the  results 
obtained  from  the  two-dimensional  linear  cacade  theory 
should  be  corrected  according  to  the  right-hand  side  term  of 
Eqn.  (4.3. 1-6).  It  is  readily  understood  that  these  right- 
hand  side  terms  are  satisfied  by  distributing  the  following 
vortices  and  sources  on  the  entire  (X,  Y)  plane 


(a)  5  =  (7  x  w)x  Y 


(b)  y  =  (V»jw)X  Y  = 


r  2 

2oj  (■£— )  sinX 
ro 

(4.3. 1-8) 

_1 ja(bp)  3jjj  3(  bp  )  3 jjj 

(bp)*'  /  ax  3 Y  3 Y  3xj  . 


By  adding  the  induced  velocities  calculated  from  c  and  y, 
the  blade  profile  shape  or  equivalently  the  camber  obtained 


in  the  conventional  two-dimensional  analysis  will  be 
corrected.  It  should  be  noted  that  the  first  term  on  the 
right-hand  side  of  Eqn.  (4.3. 1-6)  arises  from  non-zero  X, 
i.e.,  the  stream  surface  is  not  parallel  to  the  axis  of 
rotation,  whereas  the  second  group  of  terms  is  due  to  the 
non-uniform  thickness  of  stream  surface  or  tube.  Needless 
to  say,  if  X  =  0  and  bp  is  constant,  Eqn.  (4.3. 1-6)  becomes 
a  Laplace  equation  and  thus  a  two-dimensional  linear  cascade 
theory  holds. 

A  method  similar  to  the  present  one  was  developed  by  Inoue 
and  his  colleague  (e*g.,  the  paper  by  Inoue,  et  al. 
(1980)).  In  this  paper  there  exist  a  few  major  drawbacks, 
some  of  which  could  potentially  lead  to  a  substantial  error 
in  the  final  design.  First  of  all,  since  they  use  a  two- 
dimensional  linearized  cascade  theory,  the  error  becomes 
significant  for  high  solidity  and  high  stagger  angle  area, 
i.e.,  near  the  hub,  although  they  introduce  experimental 
data  in  a  later  step  of  the  analysis.  Secondly,  their  vel¬ 
ocity  triangle  used  for  determining  the  incoming  flow  angle 
to  the  blade  is  in  error  of  the  first  order  since  they  did 
not  take  into  consideration  the  effect  of  non-cylindrical 
and  variable  thickness  stream  surface.  Finally,  due  to  the 
use  of  the  linearized  cascade  theory,  they  failed  to  obtain 
the  velocity  distribution  so  that  a  boundary  layer  analysis 
and  cavitation  inception  analysis  are  not  possible. 

With  these  aspects  in  mind,  effort  has  been  made  in  the 
current  GHR  project  to  improve  the  accuracy  of  the  linear 
cascade  theory  as  well  as  to  avoid  the  singular  behavior  of 
velocity  at  the  leading  edge  of  blade,  see  Section  4.3.2  and 
Section  4.3.3,  respectively.  The  detailed  analysis  and 
numerical  procedure  for  the  vorticity  and  source  corrections 
on  the  two-dimensional  flow  are  yet  to  be  developed  and  will 
be  handled  in  the  FY-86  GHR  program. 


Loading  Correction  on  Linearized 
Cascade  Theory 

The  cascade  lift  coefficients  calculated  based  on  linearized 
theory  are  found  unsatisfactory.  Two  correction  factors  are 
proposed  in  the  present  study  to  determine  the  effective 
lift  coefficient  which  fit  better  to  the  experimental  data. 

According  to  the  linearized  cascade  theory  (Mellor,  1959), 
the  ideal  theoretical  lift  coefficient  can  be  obtained  as  a 
function  of  the  mean  angle  of  attack,  am,  camber  distribu¬ 
tion,  yc(x);  blade  thickness  distribution,  yp(x);  solidity, 
o;  and  stagger  angle,  X.  The  results  are  compared  with  the 
NACA  65-series  experimental  data  given  by  Herrig,  et_  al . 
(1951). 

The  correlation  between  the  calculated  and  measured  lift 
coefficients  is  depicted  in  Figure  4-15.  With  3i  denoting 
the  upstream  flow  angle  relative  to  the  cascade  axis,  the 
results  for  a  sample  case  of  8i  =  45’  and  a  =  1.5,  with  NACA 
65-(15)10,  are  shown  in  Figure  4-16.  These  figures  indicate 
that  the  calculated  results  obtained  from  the  linearized 
theory  are  in  general  higher  than  the  laboratory  data. 

Mellor  (1959)  did  a  similar  comparison  study  and  found  that 
the  effects  due  to  the  camber-line  slopes  at  the  blades 
extremities  may  be  suppressed  to  yield  results  closer  to  the 
laboratory  data.  Under  this  concept  of  reducing  the  camber 
effect,  he  suggested  to  take 

Cb( effective )  =  Cp( theoretical )  (4. 3. 2-1) 

when  Cp  denotes  the  single-foil  camber.  The  effective  Cp  is 
the  one  to  be  used  in  the  model  to  find  an  effective  lift 
coefficient  suitable  for  any  application. 

The  generated  lift  coefficients  calculated  following  his 
suggestion  are  compared  to  the  measured  values  and  shown  in 
Figures  4-17  and  4-18.  The  results  are  improved  when  com- 
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pared  to  Figures  4-15  and  4-16  which  are  obtained  from  the 
case  without  applying  any  modification  factor. 


After  following  Mellor’s  suggestion  of  C5  modification, 
although  the  deviation  between  the  theoretical  result  and 
laboratory  data  is  reduced,  the  obvious  difference  still 
exists,  especially  at  a  high  flow  angle.  In  order  to 
further  reduce  the  deviation,  the  present  study  employs  dual 
modification  factors  Kcb  and  Ka  such  that 


Cb(eff . )  =  KCb  Cb( theor . ) 

and 


am(eff.)  Kam  am(theor.) 


(4. 3. 2-2) 


(4. 3. 2-3) 


These  two  factors  are  determined  by  having  the  least 
deviation  between  the  calculated  lift  coefficient  and  the 
corresponding  measured  values.  The  details  are  described  in 
the  following. 

For  a  combination  of  Kcb  and  Kam,  a  set  of  as  a  function 
of  C5,  a,  81,  and  is  calculated  from  the  linearized 
model.  For  each  trial  combination  of  Kq^  and  Kam,  the  com¬ 
puted  Ci_^  is  compared  to  the  measured  values  and  a  value  of 
the  standard  deviation  is  produced.  The  set  of  Kcb  and  Kam 
which  gives  minimum  standard  deviation  is  taken  as  the 
desired  coefficients  to  produce  the  best  fitted  results. 

The  application  of  Kam  in  addition  to  K^t,  yields  further 
reduction  of  the  error.  If  the  results  are  examined  with 
the  coefficients  read  at  an  interval  of  0.05,  the  minimum 
standard  deviation  of  the  residuals  for  all  data  occurs  at 
(KCb>  Kam)  =  (0.7,  0.75).  If  the  available  data  is  divided 
into  four  groups  of  81  =  30*,  45',  60’,  and  70',  the  minima 
occur  at  (0.75,  0.8),  (0.65,  0.85),  (0.65,  0.7),  and  (0.6, 
0.75),  respectively.  If  examined  at  an  interval  of  0.01, 
these  five  sets  of  results  are  (0.69,  0.74),  (0.73,  0.83), 
(0.67,  0.79),  (0.63,  0.73),  and  (0.62,  0.71),  respectively. 


The  correlation  between  the  model  results  and  laboratory 
results  is  given  in  Figure  4-19.  The  improvements  of 
applying  the  modification  factors  can  be  found  by  comparing 
these  figures  to  Figures  4-15  and  4-17.  Without  using  the 
factors  (Figure  4-15),  the  model  overpredicts  the  results  in 
general,  with  the  higher  residual  for  higher  lift  coef¬ 
ficient  which  is  associated  with  higher  angle  of  incidence. 
After  applying  K^b  =  0.725  as  suggested  by  Mellor  (1959), 
the  standard  deviation  of  the  residuals  decreases  from  0.118 
in  Figure  4-15  to  0.060  in  Figure  4-17.  With  a  combination 
of  (KCb>  Kam)  =  C 0 . 7 ,  0.75),  the  resultant  value  fits  much 
better  to  the  measured  value  (Figure  4-19)  and  the  standard 
deviation  becomes  only  0.026. 

With  only  j-,  =  0.725,  a  sample  comparison  of  the  calculated 
and  measured  lift  coefficient  (Figure  4-18)  shows  an 
improvement  from  that  in  Figure  4-16.  Yet  the  results  are 
still  too  high  and  are  worse  at  higher  angle  of  attack. 
With  (Kcb»  Kam)  =  (0.7,  0.75),  the  result  is  much  better,  as 
shown  in  Figure  4-20.  The  residuals  are  now  more  or  less 
evenly  distributed  over  the  angle  of  attack  (and  over  the 
magnitude  of  lift  coefficient). 

This  study  of  loading  correction  on  linearized  cascade 
theory  is  conducted  to  provide  realistic  loading  data.  Both 
the  single-foil  camber,  Cb,  and  angle  of  attack,  am,  are 
multiplied  by  discount  factors  before  they  are  used  in  the 
theoretical  model  to  calculate  the  cascade  lift  coefficient. 
The  discount  factors  are  determined  based  on  a  set  of 
laboratory  data.  This  engineering  approach  of  obtaining  the 
desired  lift  coefficient  is  deemed  appropriate  to  provide 
data  to  be  used  in  a  propulsor  blade  design  procedure. 

In  the  linearized  cascade  theory  (Mellor,  1959),  it  is 
assumed  that  the  difference  between  the  induced  velocity  on 
the  blade  surface  and  on  the  camber  line  is  negligible. 
This  assumption  is  not  valid  under  certain  circumstances. 
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The  errors  are  proportional  to  the  square  of  camber  and 
square  of  relative  thickness  as  indicated  by  Mellor  (1959). 
Thus  it  is  planned  to  take  a  theoretical  approach  in  the 
future  study,  instead  of  the  present  engineering  approach, 
to  correct  the  errors  due  to  assumptions  adopted  in  the 
linearized  theory.  The  singularities  are  to  be  distributed 
on  the  camber  line  or  even  on  the  blade  surfaces  instead  of 
on  the  chord  line.  In  this  way,  the  results  will  be  more 
accurate  by  paying  for  extra  complicity  in  the  solution  pro¬ 
cedure. 


4.3.3  Leading  Edge  Correction  on  Linearized 
Cascade  Theory 

The  second  disadvantage  of  using  the  linearized  cascade 
theory  stems  from  the  fact  that  a  singularity  in  terms  of 
flow  velocity  exists  at  the  blade  leading  edge.  This  singu¬ 
larity  is  integrable  so  that  it  will  present  no  problem  of 
calculating  the  force  on  the  blade.  It  will,  however,  cause 
inconvenience  whenever  a  detailed  velocity  distribution  is 
required,  e.g.,  for  calculation  of  viscous  boundary  layer  or 
cavitation  inception  on  the  blade. 

The  objective  of  the  study  in  this  Section  is  therefore  to 
correct  such  problem  and  then  to  provide  an  accurate  veloc¬ 
ity  distribution.  The  method  employed  here  is  a  singular 
perturbation  method  widely  used  in  fluid  flow  problems,  see 
e.g.,  a  textbook  by  Van  Dyke  (1975),  a  paper  by  Furuya  and 
Acosta  (1973).  The  method  is  particularly  useful  for  a  case 
in  which  a  solution  is  regular  everywhere  except  for  a 
localized  singularity. 

The  velocity  obtained  from  a  linearized  cascade  theory  (see 
Mellor  (1959))  is  given 

oo 

=  1  +  2ttAq-  rQ0  +  2ttA  r01  +  2  X  A  n  sin  ne  +  p(  e(x) ) 

n=2 

(4.3.3-la) 
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where  +  corresponds  to  the  upper  and  lower  side  of  the 
blade, 

^  r  2fl0  ^sVna  *  4  nf>n  sln  n8> 

qm  =  geometric  mean  velocity  of  upstream  and 
downstream  velocities,  and  q2, 

cos  9  =  1-2  x/c, 

cos  eQ  =  1-2  xq/c,  where  9Q,  xQ  are  used  as  dummy 
variables , 

c  =  chord  length, 

yc  =  Cb>fVc^  camber  Function ) , 
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X  =  stagger  angle  of  cascade, 
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and  , 
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As  clearly  seen  from  I^g  term,  the  velocity  becomes  singular 
at  the  leading  edge,  i.e.,  as  x  +  0  (i.e.,  e  -*■  0).  The  lift 
coefficient  can  be  calculated  based  on  this  velocity 


Pp-P, 


Ci  *  T 


2  H  m 


=  2ttA 0  +  2uA  r  (4. 3. 3-2) 

It  is  noted  that,  although  the  velocity  has  a  singularity  at 
the  leading  edge,  the  velocity  squared  is  integrable  so  that 
the  force  can  be  conveniently  calculated.  It  should  also  be 
mentioned  that  only  the  first  two  circulation  terms,  i.e., 

-Hdo ,  ,  remain  for  the  lift  calculation.  By  expanding  the 
leading  edge  area  of  y  =  +  e/  x  with  x  =  e2X  and  y  =  e2Y 
where  the  leading  edge  radius  R  =  e2/2 ,  the  flow  velocity 
becomes  (see  the  paper  of  Furuya  (1983)), 


qi  =  Ui  /x  +  l/4 


(4. 3. 3-3) 


where  X=a2  denotes  the  X-coordinate  of  stagnation  point  and 
Uj_  is  the  flow  velocity  at  upstream  infinity  in  the  inner 
region.  The  inner  region  is  represented  by  (X,Y)  which  is 
stretched  by  the  scale  of  e2.  Also,  in  the  present  method, 
it  is  assumed  that  the  round  leading  edge  profile  of  any 
blade  can  be  accurately  approximated  by  a  parabola. 
Rewriting  Eqn.  (4. 3. 3-3)  in  terms  of  x  and  expanding  it  as 

X  ■*  <», 


qi  =  Ui 
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=  U,  (l  +  -£2_ 
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+  0  ( e2 ) 
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(4. 3. 3-4) 


Expansion  of  the  outer  solution  in  Eqn.  (4.3.3-la)  as  x+0 
becomes,  after  some  algebra, 


q~  =  1  - 

“  m 


+  0(e). 


(4. 3. 3-5) 


Matching  the  inner  solution  (4. 3. 3-4)  with  the  outer  solu¬ 
tion  (4. 3. 3-5)  yields 


Ui  =  dm 


ea  =  A, 


(4. 3. 3-6) 


The  uniformly  valid  solution  (see  the  book  of  Van  Dyke 
(1975))  can  then  be  constructed 


qUqlf°r"'  *  1  i  2wA0r00c»)  i  2”A1  r01(9) 
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+  2  22  fln  sin  n®  + 


x+e2/4 


(4. 3. 3-7) 


where  F(e(x))  is  defined  in  Eqn.  (4.3.3-lb).  It  should  be 
noted  that  the  uniformly  valid  solution  no  longer  has  a 
singularity  since  the  second  term  together  with  the  last 
term  in  Eqn.  (4. 3. 3-7)  now  has  a  finite  value  as  x-*-0. 

The  above  results  were  used  to  compare  the  analytically  pre¬ 
dicted  velocity  profile  with  experimental  data.  Four  dif¬ 
ferent  such  results  using  NACA  65-series  airfoil  are  shown 
in  Figure  4-21.  Not  only  the  singular  behavior  which  would 
exist  with  the  linearized  cascade  theory  has  disappeared  but 
also  the  agreement  is  good  enough  to  be  applicable  for 
boundary  layer  and  cavitation  inception  analyses. 
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CONCLUSIONS 


Through  the  FV-85  GHR  program,  the  following  major  subjects 
have  been  investigated: 

o  It  has  been  concluded  that  among  three  possible  can¬ 
didate  design  theories  which  can  incorporate  three- 
dimensionality  into  pumpjet  design  procedure,  the 
blade-through  flow  with  blade-to-blade  flow  method 
has  been  selected  as  the  appropriate  design  method 
for  the  present  study. 

o  The  selected  design  method  is  of  an  iterative  type, 
which  includes  the  following  sub-theories: 

(a)  the  determination  of  the  shroud  intake  diameter, 

(b)  the  determination  of  the  streamline  by  a  stream¬ 
line  curvature  method  (SCM), 

(c)  mapping  of  the  stream  surface  calculated  through 
the  SCM  onto  a  plane  hence  producing  a  Poisson 
equation  due  to  the  cone  shaped  stream  surface, 

(d)  modifying  the  blade  profile  shape  obtained  in  the 
potential  theory  through  a  method  of  correcting 
the  effect  of  the  Poisson  equation  on  the  poten¬ 
tial  theory  results, 

(e)  improving  the  prediction  accuracy  of  the  two- 
dimensional  cascade  theory  by  using  a  pair  of 
correction  factors  such  that  the  loadings  calcu¬ 
lated  through  the  linearized  cascade  theory  fit 
well  with  the  experimental  data. 

(f)  removing  the  singularity  at  the  leading  edge  of 
the  two-dimensional  cascade  theory  by  applying 
the  singular  perturbation  method  so  that  a  useful 
velocity  profile  can  be  obtained  to  be  later  used 
for  flow  separation  and  cavitation  analyses. 

(g)  checking  for  the  possibility  of  flow  separation, 
and 

(h)  checking  for  the  chance  of  cavitation  inception. 


o  Among  the  above  sub-theories  (a)  -  (h),  the  theory 

development  of  (a),  (b),  (e),  and  (f)  has  been 

accomplished  during  FY-85. 

For  the  FY-86  GHR  program,  it  is  planned  to  perform  the 
following : 


o  develop  the  detailed  theoretical  formula  for  the 
diagonal  flow,  (c)  and  (d)  above, 

41  o  develop  a  theoretical  approach  to  improve  the  result 

of  loading  calculated  through  a  linearized  cascade 
theory,  e.g.,  having  singularities  distributed  on  the 
camber  instead  of  the  chord  line,  an  additional  work 
for  (e),  and 

•  o  develop  a  foundation  for  a  unified  pumpjet  design 

theory  by  incorporating  all  the  sub-theories  devel¬ 
oped  above. 
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APPENDIX 


CALCULATION  OF  HEAD  LOSS  COEFFICIENT  Kx 
(Between  Reference  Point  (T)  and  Rotor  Inlet  (2) ) 

The  energy  loss  between  (T)  and  (2)  can  be  given 
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Defining  Pt  =  Ps  +  2pV  (total  pressure), 
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(total  pressure  coefficient), 


(static  pressure  coefficient) 
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Since  Cpt^,  CpS2  and  V2  are  given  as  a  function  of  r  in  the 

report  of  Treaster  (  1977)  (see  Figures  4  and  5),  we  can 
2 

determine  AH/y  and  thus  as  follows 
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AH 
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v„/2g 


(7) 


2g 


v, 

v~ 


2  »  1 

where  AH/V  /  ,  and  —  (Eqn.  (24))  are  both  calculated. 
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Qualitative  Comparisons  for  the 
Candidate  Methods  as  a  Three-Dimensional 
Pumpjet  Design  Method 
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FIGURE  1-4  Typical  pump jet  rotor  blade  con¬ 
figuration,  (a)  top  view  and 
(b)  upstream  view 
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FIGURE  1-5  Typical  Cp  and  Gp  curves  for  the 
designed  rotor  blade  profiles  and 
deformed  rotor  blade  profiles, 
showing  about  7.3%  increase  of  Cp 
(i.e.,  ACd  =  .0083)  where  the  design 
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Stream  Surface  by 
Streamline  Curvature 
Method  with  Blade 
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FIGURE  4-1  Flow  Chart  of  the  Selected 
Pumpjet  Design  Method 
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FIGURE  4-2  Schematic  diagram  of  pumpjet  flow 


Cm  can  be  calcu¬ 
lated  from  (4.1.1-14) 
and  (4.1.1-15)  or 
(4.1.1-10) 


ei,  07  given,  Vi 
can  be  calculated 
from  (4. 1.1-5) 


Ki  can  be  calcu¬ 
lated  from  the 
formula  in  Appen¬ 
dix;  Vi  can  be  cal¬ 
culated  from  (4. 1.2-5) 


fir  should  be  given, 
usually  riR  ~  89% 


FIGURE  4-3  Flow  chart  for  calculation 
of  the  pumpjet  efficiency 


FIGURE  4-5  Meridional  \/elocity  distributions  at 
stations  @  and  (3)  of  Figure  4-2 
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FIGURE  4-9  Maximum  propulsive  efficiency  as  a  function 
of  the  incoming  flow  velocity  amplification 
factor  with  0-,  as  a  parameter 


FIGURE  4-10  Velocity  diagram 
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FIGURE  4-11  A  schematic  flow  diagram  used 
for  numerical  computations  on 
Streamline  Curvature  Method 
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FIGURE  4-15  Comparison  of  lift  coefficient  data, 
without  modification  factor 
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FIGURE  4-16  Comparison  of  theoretical  (dashed 
line)  and  measured  (solid  line) 
lift  coefficient  data  for 
C  NACA  6  5- ( 15)  10 ,  6,  =  45°, 

a  =  1.5  without  modification 
factor.  (The  arrow  indicates 
the  design  angle  of  attack.) 
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©  FIGURE  4-17  Comparison  of  lift  coefficient  data, 

with  a  modification  factor  of  0.725 
applied  to 


© 


68 


0 . 3 


1 


12 


16 


20 


J _ L 


24 


28 


32 


FIGURE  4-18 


Comparison  of  theoretical  (dashed 
line)  and  measured  (solid  line)  lift 
coefficient  data  for  NACA  65- (15) 10, 
Bi  =  45°,  o  =  1.5,  with  Kcb  =  0.725 
and  K  =  1.  (The  arrow  indicates 
the  diSign  angle  of  attack.) 


69 


CALCULATED  LIFT  COEFFICIENT 


MEASURED  LIFT  COEFFICIENT 


6 


FIGURE  4-19  Comparison  of  lift  coefficient  data, 
with  (Kcb,  Kam)  =  (0.7,  0.75) 
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FIGURE  4-20  Comparison  of  theoretical  (dashed  line) 
and  measured  (solid  line)  lift  coef¬ 
ficient  data  for  NACA  65- (15) 10,  3^  = 
45°,  a  =  1.5,  with  Kcb  =0.7  and 
Kccm  =  °’75*  (The  arrow  indicates 
the  design  angle  of  attack.) 
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FIGURE  4-21  Comparison  of  velocity  distribution  between 

the  singular  perturbation  method  ( - )  and 

test  results  (0:  upper  surface,  ■:  lower 
surface)  of  Herrig,  et  al.  (1951)  for 
a)  NACA  65- (12) 10,  0,  =  45°,  solidity  = 

1.0  and  a1  =  12.1°  and  b)  NACA  65-(12)10, 

8^  =  45° , isolidity  =  0.5  and  =  7.3° 
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FIGURE  4-21  (cont’d)  c)  NACA  65-410,  Bt 

=  45°, 

•solidity  =  1.0  and 

oi  =  5°, 

and  d)  NACA  65-810, 

B,  =  45 

solidity  =  1.0  and 

aL=  9.7 
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